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1.  SYNOPSIS  OF  CONTRACT  ACTIVITIES 

During  the  past  seven  years,  we  have  made  substantive  contributions  toward  the 
goals  set  forth  in  our  original  grant  proposal.  In  the  research  area,  we  have  succeeded  in 
developing  the  basic  understanding  of  several  mesoscale  space  plasma  phenomena  of 
considerable  importance.  These  include  the  study  of  the  origin  of  high-latitude 
magnetospheric-ionospheric  turbulence,  the  formation  of  ion  conics  in  the  auroral  and 
cusp  regions,  the  generation  of  counterstreaming  electrons  along  auroral  field  lines,  the 
process  leading  to  the  anomalous  heat  transport  in  the  polar  wind  due  to  photoelectrons, 
the  development  of  black  auroral  curls  in  the  return  current  region,  the  apparent  low¬ 
dimensional  behavior  of  self-organized  criticality  and  global  geomagnetic  activity,  and 
other  related  nonlinear  space  plasma  processes.  A  number  of  these  theoretical  findings 
have  recently  been  confirmed  by  experimental  data  collected  by  high-time  resolution 
instruments  on  board  the  Freja,  Akebono,  POLAR  and  FAST  satellites  as  well  as  the 
TOPAZ,  SCIFER  and  AMICIST  rocket  series.  As  a  gauge  of  the  seminal  nature  and 
excellence  of  our  contributions  and  their  popular  acceptance  by  our  peers,  we  note  that  we 
were  invited  to  deliver  a  total  of  over  50  invited  and  review  lectures  at  various  national 
and  international  conferences  and  major  prestigious  institutions.  We  were  also  asked  to 
chair  and/or  organize  specialty  sessions  on  particle  acceleration,  polar  wind  and  multiscale 
phenomena  at  a  number  of  conferences  and  workshops.  During  the  contract  period,  we 
published  a  total  of  40  scientific  papers  and  five  proceedings  volumes. 

We  have  interacted  with  a  number  of  research  organizations  including  the 
Air  Force  Research  Laboratory,  the  Naval  Research  Laboratory,  Cornell  University,  the 
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International  Space  Science  Institute,  the  Universities  of  California  at  Berkeley,  Los 
Angeles  and  Irvine,  the  University  of  Maryland,  the  University  of  New  Hampshire,  Boston 
College,  the  Max-Planck  Institutes  for  Extraterrestrial  Physics  and  Aeronomy,  the 
Lockheed  Palo  Alto  Research  Laboratory,  the  National  Research  Council  of  Canada,  the 
University  of  Iowa,  the  Swedish  Space  Institute  and  the  Space  Research  Institute  in 
Moscow.  Visits  to  these  institutions  and  by  scientists  from  these  institutions  have 
provided  the  necessary  stimulus  to  keep  our  research  program  vibrant,  up-to-date,  and  at 
the  same  time  constantly  in  touch  with  practical  motivations. 

We  organized  several  combined  Symposia/Workshops  on  the  "Physics  of  Space 
Plasmas"  covering  the  topics  of  "Chaos,  Stochasticity,  and  Strong  Turbulence."  and 
"Multiscale  Phenomena  in  Space  Plasmas."  In  addition.  Dr.  Tom  Chang  organized  a 
Memorial  Symposium  for  the  Nobel  Laureate,  Professor.  S.  Chandrasekhar  of  the 
University  of  Chicago.  This  activity  was  held  in  conjunction  with  the  Annual  Meeting  of 
the  Society  of  Engineering  Science  in  October  1996  at  the  Campus  of  the  Arizona  State 
University  in  Tempe.  This  event  attracted  over  300  scientists  and  engineers  world-wide. 

The  content  of  this  report  is  organized  as  follows.  In  Sections  2  -  4,  we  discuss  in 
detail  the  three  most  noteworthy  contributions  of  our  research  activities:  (1)  ion 
energization  and  outflow,  (2)  self-organized  criticality  and  multiscale  intermittent 
turbulence  in  the  plasma  sheet,  and  (3)  a  new  theory  of  the  nonclassical  polar  wind.  In 
Section  5,  a  complete  listing  of  the  scientific  publications  of  the  MIT  group  during  the 
contract  period  (including  the  titles  of  five  edited  conference  proceedings)  is  provided. 
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Ion  Energization  and  Outflow: 

Some  Important  Mechanisms 

Mats  Andre  and  Patrik  Norqvist 

Swedish  Institute  of  Space  Physics,  Department  of  Theoretical  Space  Physics,  Umea  University,  901  87  Umea,  Sweden 


Tom  Chang 

Center  for  Space  Research,  Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139,  USA 


Abstract.  More  than  one  mechanism  is  causing  ion 
energization  perpendicular  to  the  geomagnetic  field  in 
the  terrestrial  auroral  region.  Ions  energized  in  this 
way  may  leave  the  ionosphere,  and  are  an  important 
source  of  magnetospheric  plasma.  Observations  by  the 
Freja  satellite  obtained  during  20  months  at  altitudes 
around  1700  km  are  used  to  study  the  most  common 
0+  energization  mechanisms.  The  most  common  and 
most  important  ion  energization  is  due  to  broadband 
low-frequency  waves,  causing  resonant  energization  by 
interaction  at  frequencies  around  the  ion  gyrofrequency. 
This  mechanisms  is  compared  with  other  types  of  ion 
energization.  The  dependence  of  the  different  energiza¬ 
tion  types  on  magnetic  local  time  and  magnetic  latitude 
is  presented.  We  also  speculate  on  the  origin  of  the 
broadband  waves  causing  most  of  the  ion  energization. 


1.  Introduction 

An  intriguing  problem  of  magnetospheric  physics  is 
the  energization  of  ionospheric  ions  and  their  subse¬ 
quent  escape  from  low  altitudes.  Much  of  the  ion  out¬ 
flow  is  caused  by  ion  heating  in  the  direction  perpendic¬ 
ular  to  the  geomagnetic  field.  The  heated  ions  may  then 
move  up  the  field  lines  of  the  inhomogeneous  terres¬ 
trial  magnetic  field  and  form  so-called  conics  in  velocity 
space.  This  ion  energization  is  often  associated  with  the 
auroral  region  and  occurs  at  all  local  times.  Several  ion 
energization  mechanisms  have  been  suggested,  but  the 
relative  importance  of  the  different  mechanisms  has  not 
previously  been  established.  Using  observations  of  0+ 
ions  obtained  by  the  Freja  satellite,  we  present  a  sta¬ 
tistical  study  of  the  most  common  0+  ion  energization 
mechanisms. 

Observations  of  perpendicular  ion  heating  have  been 
obtained  during  several  years  from  rockets  and  satel¬ 
lites  at  altitudes  from  a  few  hundred  kilometer  and  at 
least  up  to  several  Earth  radii.  The  ion  distributions 
may  have  average  energies  from  a  few  eV  up  to  at  least 
a  few  keV,  i.e.,  thousands  of  times  higher  than  the  mean 


ion  energies  in  the  lower  ionosphere.  Most  of  the  ener¬ 
gization  takes  place  above  the  collisional  ionosphere  and 
thus  depends  on  interaction  between  the  charged  ions 
and  some  electric  fields.  The  importance  of  ionospheric 
ion  outflow  as  a  source  of  magnetospheric  plasma  is 
well  established  by  statistical  studies  using  polar  or¬ 
biting  satellites,  as  discussed  in  reviews  by  Moore  and 
Delcourt  [1995]  and  Yau  and  Andre ,  [1997].  Several  ion 
energization  mechanisms  causing  this  outflow  have  been 
suggested  (see  reviews  by,  e.g.,  Klumpar,  1986;  Lysak, 
1986,  Andre  and  Chang ,  1993;  Andre ,  1997;  Andre  and 
You,  1997;  Chang  and  Andre)  1998) 

Recent  observations  indicate  that  more  than  one  ion 
energization  mechanism  is  operating  in  the  terrestrial 
auroral  region,  but  also  suggest  that  one  mechanism  is 
the  most  important.  Broadband  low-frequency  electric 
wave  fields  are  commonly  observed  together  with  per¬ 
pendicular  ion  energization.  These  waves  cover  frequen¬ 
cies  from  below  one  Hz  up  to  several  hundred  Hz,  thus 
including  the  gyrofrequencies  of  the  major  ion  species 
at  least  for  altitudes  from  about  1000  km  up  to  sev¬ 
eral  Earth  radii.  Several  wave  and  particle  observations 
by  satellites  [e.g.,  Norqvist  et  al}  1996;  Carlson  et  al. , 
1998]  and  by  rockets  [e.g.,  Lynch  et  al. ,  1996;  Kintner 
et  al.,  1996]  indicate  a  clear  correlation  between  these 
waves  and  ion  energization.  Resonant  energization  by 
the  observed  waves  at  frequencies  near  the  ion  gyrofre¬ 
quency  usually  can  heat  the  ions  to  the  observed  ener¬ 
gies  [Chang  et  a/.,  1986;  Retterer  et  al .,  1987;  Jasperse , 
1998].  In  the  following  we  give  an  overview  of  some  re¬ 
cent  statistical  studies  of  observations  obtained  by  the 
Freja  satellite.  We  conclude  that  resonant  heating  by 
broadband  low-frequency  waves  is  the  dominating  ion 
energization  mechanism. 

2.  An  Example  of  Ion  Energization  Observed  by 
the  Freja  Satellite 

The  Freja  satellite  was  launched  October  6,  1992, 
into  an  orbit  with  63°  inclination,  an  apogee  in  the 
northern  hemisphere  of  1750  km,  and  a  perigee  of  600 
km.  The  Freja  instruments  include  the  TICS  ion  mass 
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spectrometer,  the  TESP  electron  spectrometer,  experi¬ 
ments  for  observing  low  frequency  electric  and  magnetic 
fields,  and  experiments  for  electric  and  magnetic  wave 
fields  and  the  plasma  density.  The  Freja  instruments 
are  described  in  Space  Science  Reviews,  70 ,  405-602, 
1994.  Observations  by  Freja  have  been  used  to  de¬ 
fine  a  few  different  types  of  ion  energization.  Below 
we  present  an  example  of  an  ion  energization  event. 

An  intense  morningside  ion  heating  event  occurred 
on  December  3,  1993,  on  orbit  5591.  The  corrected  ge¬ 
omagnetic  latitude  (CGL)  was  about  72°,  the  magnetic 
local  time  (MLT)  was  around  0600,  and  the  Kp  index 
was  5.  Figure  1  shows  an  overview  of  the  observations 
(see  also  Eliasson ,  [1996]  and  Andre  et  a/., [1998]).  Note 
that  Freja  is  moving  essentially  tangential  to  the  auro¬ 
ral  oval  (at  nearly  constant  CGL)  during  the  event. 

Oxygen  ion  count  rates  from  the  TICS  ion  mass 
spectrometer  are  shown  in  panel  1  of  Figure  1.  The 
corresponding  pitch  angles  are  displayed  in  panel  3. 
Heated  ions  can  be  seen  as  vertical  stripes  in  panel  1 
between  0842:10  and  0842:30,  and  also  between  0843:10 
and  0844:55  UT.  It  is  the  heating  of  these  0+  ions  we 
consider.  Panel  2  shows  data  from  the  TICS  proton 
channel,  which  is  not  used  in  our  statistical  studies. 

Panel  4  displays  downgoing  electrons  observed  by 
the  TESP  detector.  Note  that  the  downgoing  keV  (so- 
called  inverted  V)  electrons  are  anticorrelated  with  the 
ion  energization.  We  have  found  that  during  several 
events,  on  a  timescale  of  minutes  (several  hundreds  of 
kilometers  along  the  satellite  trajectory)  the  ion  heat¬ 
ing  and  keV  auroral  electrons  are  associated,  but  on  a 
timescale  of  seconds  (tens  of  kilometers)  the  ion  ener¬ 
gization  and  the  auroral  electrons  are  anticorrelated. 

Panel  6  shows  a  time  series  of  the  perpendicular  low- 
frequency  (below  tens  of  Hz)  electric  field  in  the  satel¬ 
lite  spin  plane,  while  panel  8  displays  a  time  series  of 
the  low-frequency  magnetic  field  fluctuations.  The  am¬ 
plitude  of  these  electric  and  magnetic  fluctuations  are 
correlated  with  the  0+  ion  energization. 

Panel  5  displays  an  electric  field  spectrum  up  to  16 
kHz.  There  is  a  good  correlation  between  waves  with 
frequencies  from  below  one  hertz  up  to  several  times  the 
oxygen  gyrofrequency  /0+  at  25  Hz  and  the  ion  heating. 
There  are  also  some  waves  in  the  lower  hybrid  frequency 
range  (1  to  2  kHz),  but  there  is  no  detailed  correlation 
with  the  ion  heating.  Panel  7  shows  a  spectrum  of  the 
magnetic  field  below  64  Hz.  The  magnetic  spectral  den¬ 
sity  is  also  correlated  with  the  ion  heating.  However, 
the  magnetic  spectrum  decreases  more  rapidly  with  in¬ 
creasing  frequency  than  the  electric  spectrum.  Panel  9 
displays  the  density,  obtained  from  a  Langmuir  probe. 

To  summarize,  a  clear  correlation  between  the  trans¬ 
versely  heated  oxygen  ions  and  waves  with  frequencies 
from  below  one  hertz  up  to  several  times  /0+  is  ob¬ 
served  for  this  event.  It  can  be  shown  that  generally 
the  observed  wave  amplitudes  are  high  enough  to  heat 
the  0+  ions  to  the  observed  energies  for  such  events 
[Andre  et  al. ,  1998]. 


3.  Different  Types  of  Ion  Energization 

By  considering  in  detail  data  sets  obtained  by  Freja 
during  20  events  at  various  local  times,  a  few  different 
types  of  ion  heating  events  can  be  defined.  An  example 
of  a  data  set  is  shown  in  Figure  1. 

Type  1  ion  heating  is  by  Norqvist  et  ai  [1998]  defined 
to  be  associated  with  broadband  low-frequency  waves, 
similar  to  the  event  in  Figure  1.  An  increase  of  the 
field-aligned  current  is  often  associated  with  this  type 
of  ion  heating.  The  heating  may  also  be  associated 
with  suprathermal  electron  bursts,  i.e.,  field  aligned 
electrons  covering  a  broad  energy  range  from  tens  of 
eV  up  to  several  hundred  eV. 

Type  2  ion  heating  is  also  associated  with  broadband 
low-frequency  waves,  although  the  ratio  of  the  electric 
and  magnetic  field  is  typically  smaller  than  for  type  1 
events.  These  events  are  associated  with  precipitating 
protons  with  keV  energies,  and  electrons  with  mean  en¬ 
ergies  of  about  100  eV,  typical  of  the  cusp/cleft  region. 
Similar  to  type  2  events,  type  1  events  are  often  associ¬ 
ated  with  field-aligned  currents,  and  with  suprathermal 
electron  bursts. 

Type  3  ion  heating  is  associated  with  downgoing  keV 
auroral  electrons.  Type  3(LH)  events  are  associated 
with  waves  near  the  lower  hybrid  frequency,  while  type 
3(EMIC)  are  associated  with  electromagnetic  ion  cy¬ 
clotron  waves  at  about  half  the  proton  gyrofrequency. 

4.  A  Statistical  Study  of  Ion  Energization  Mech¬ 
anisms 

Using  Freja  observations  obtained  in  the  northern 
hemisphere  from  late  January  1993  to  the  middle  of 
September  1994,  a  statsitical  study  of  ion  energization 
mechanisms  has  been  performed  [Norqvist  et  al .,  1998]. 
This  period  coincides  with  the  declining  phase  of  solar 
cycle  22.  The  corrected  geomagnetic  latitude  (CGL)  of 
the  events  is  from  about  50°  to  75°,  the  upper  limit 
being  set  by  the  63°  degree  inclination  of  the  orbit.  Es¬ 
sentially  all  observations  are  obtained  between  1400  and 
1750  km,  with  a  majority  of  the  observations  ocurring 
around  1700  km. 

Since  some  of  the  upward  flowing  ions  may  come 
back  to  low  altitudes,  we  use  the  term  upflow,  rather 
than  outflow,  in  the  following.  This  upflow  is  calcu¬ 
lated  as  upward  ion  flux  minus  downward  ion  flux,  and 
the  flux  is  corrected  for  the  satellite  motion  parallel  to 
the  geomagnetic  field.  The  upflow  is  normalized  with 
respect  to  the  total  number  of  satellite  spins  in  the  ap¬ 
propriate  bin,  using  zero  flux  for  spins  not  regarded  as 
an  event.  Finally  the  upflow  is  projected  to  a  reference 
altitude  of  1000  km. 

Useful  data  were  obtained  from  about  700  000  six 
second  spacecraft  spins,  of  which  about  500  000  oc¬ 
curred  at  invariant  latitudes  above  60°.  In  the  following 
statistics,  ion  energization  of  types  3(LH)  and  3(EMIC) 
are  treated  together  as  type  3.  The  relative  importance 
of  lower  hybrid  and  EMIC  waves  is  discussed  in  section 
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5.  A  total  of  5211  spins  with  heated  0+  ions  includes 
3358  (64.4  %)  of  type  1,  1213  (23.3  %)  of  type  2,  and 
640  (12.3  %)  of  type  3.  As  demonstrated  later,  type  2 
is  associated  with  higher  upflow  and  also  dominates  the 
total  upflow. 


4.1.  Magnetic  Local  Time  Dependence 

The  occurrence  frequency  of  each  heating  type  as  a 
function  of  magnetic  local  time  for  CGL>60°  is  shown 
in  Figure  2a.  Type  1  is  the  most  common  heating  type 
on  the  nightside.  These  events  are  distributed  around 
midnight  with  a  slight  predominance  on  the  morning- 
side.  Type  2  is  the  most  common  energization  type 
around  noon.  These  events  are  significantly  more  com¬ 
mon  in  the  prenoon  sector,  as  compared  to  the  post- 
noon  sector.  The  type  3  events  are  mainly  found  at  the 
nightside.  These  events  are  clearly  more  common  on 
the  evening  side  than  on  the  morningside. 

The  average  0+  upflow  (upward  number  flux)  as  a 
function  of  MLT  for  all  three  heating  types,  is  shown 
in  Figure  2b.  Type  2  events  in  the  prenoon  sector  give 
the  highest  upward  0+  flux.  One  reason  that  type  2  on 
the  dayside  dominate  the  upflow,  although  type  1  on 
the  nightside  has  higher  occurrence  frequency,  is  that 
the  plasma  density  is  higher  on  the  dayside.  The  av¬ 
erage  upflow  during  a  spin  regarded  as  an  event  can 
be  obtained  by  dividing  the  upflow  in  Figure  2b  with 
the  occurrence  frequency  in  Figure  2a.  Thus,  the  aver¬ 
age  0+  upflow  during  a  prenoon  event  is  a  few  times 
1012/(m2s),  but  large  upflows  can  be  at  least  ten  times 
higher  [Andre  et  al ,  1998].  A  maximum  particle  upflow 
in  the  prenoon  sector  agrees  with  many  previous  stud¬ 
ies,  [see  e.g.,  Lockwood  et  a/.,  1985].  About  90  %  of  the 
events  in  Figure  3  are  of  type  1  or  2,  and  about  95  %  of 
the  0+  upflow  is  related  to  these  events.  Since  the  type 
1  and  2  are  associated  with  broadband  low-frequency 
waves,  we  find  that  most  of  the  ion  upflow  is  caused  by 
these  waves. 


Figure  2.  a)  Occurrence  frequency  for  the  three  ion  heating 
types  versus  MLT.  b)  The  corresponding  average  upward 
0+  flux  (averaging  include  all  samples  with  zero  upward 
flux)  projected  to  an  altitude  of  1000  km.  Type  1  dom¬ 
inates  the  nightside,  while  type  2  is  most  common  on  the 
dayside.  Both  types  are  caused  by  broadband  low-frequency 
waves,  and  together  completely  dominate  the  total  0+  up¬ 
flow.  (From  Norqvist  et  aL  [1998].) 


The  occurrence  frequency  for  oxygen  heating  events 
of  type  My  /m  ,  we  define  as 


/m  (*,  j) 


N(i,j) 


(1) 


where  tim  is  the  number  of  samples  (satellite  spins)  of 
type  M  events  and  N  is  the  total  number  of  samples 
in  a  specific  bin.  The  binning  parameters  (MLT,  CGL) 
are  denoted  i  and  j.  At  the  altitude  and  latitude  of 
our  study  the  difference  between  CGL  and  INVL  (used 
in  the  invariant  geomagnetic  dipole  coordinate  system) 
is  small,  and  can  be  neglected.  In  the  data  presen¬ 
tation,  summation  of  one  of  the  indices  is  performed. 
This  summation  is  made  over  all  values  of  the  binning 
parameter  if  nothing  else  is  stated. 


4.2.  Magnetic  Latitude  Dependence 

The  occurrence  frequency  of  each  heating  type  as 
a  function  of  corrected  magnetic  latitude  is  shown  in 
Figure  3a.  Type  1  has  a  peak  in  occurrence  frequency 
around  70°,  and  type  3  peaks  at  slightly  lower  latitudes. 
In  contrast,  the  occurrence  frequency  of  type  2  increases 
up  to  75°,  the  highest  latitude  included  in  the  study, 
indicating  a  maximum  at  even  higher  latitudes.  The 
upflow  in  Figure  3b  shows  a  dependence  on  latitude 
similar  to  that  of  occurrence  frequency  for  all  types. 
The  peak  in  upflow  associated  with  type  1  ion  heating 
at  62°  is  caused  mainly  by  data  from  a  single  orbit 
with  many  spins  of  high  density  conics,  and  is  believed 
not  to  be  statistically  significant.  As  in  Figure  2b  we 
find  that  the  highest  upflows  are  caused  by  type  2  ion 
heating.  An  upflow  increasing  with  latitude  in  Figure 
3b  is  consistent  with  results  from  the  DE-1  satellite 
which  show  a  peak  in  upward  0+  flux  near  noon,  where 
type  2  events  dominate,  around  78°  (Figure  8  of  Yau  et 
al  [1985]).  However,  the  DE-1  peak  flux  near  midnight, 
where  type  1  and  type  3  events  dominate,  is  just  below 
70°,  again  consistent  with  Figure  3b. 
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Figure  3.  a)  Occurrence  frequency,  and  b)  upfiow  associ¬ 
ated  with  three  ion  heating  types,  as  a  function  of  Corrected 
Geomagnetic  Latitude.  (From  Norqvist  et  al.  [1998].) 


4.3.  Associated  Electron  Energy 

The  energy  flux  of  field-aligned  downgoing  electrons 
in  and  near  regions  of  ion  energization  were  studied  by 
Andre  et  al.  [1998].  During  type  1  events  (e.g.,  Figure 
1)  it  is  clear  that  high  electron  energies  are  found  near, 
but  outside,  the  intense  ion  energization.  Figure  4a 
shows  an  anticorrelation  of  electron  energy  flux  and  ion 
energy.  Recent  observations  by  the  FAST  satellite  show 
that  very  field-aligned  electron  beams  reaching  energies 
of  hundreds  of  eV  may  be  common  in  type  1  ion  heating 
regions  [Carlson  et  a/.,  1998].  Some  of  these  very  field- 
aligned  particles  may  sometimes  be  missed  by  the  Freja 
detectors,  but  their  mean  energy  is  typically  lower  than 
that  of  auroral  electrons  detected  nearby.  During  type 
3  events,  high  electron  energy  fluxes  are  found  to  be 
directly  associated  with  the  ion  energization,  Figure  4b. 
The  relation  between  the  electron  fluxes  and  the  waves 
causing  the  ion  energization  is  discussed  below. 


5.  Waves  Associated  with  Ion  Energization 

We  first  consider  the  waves  associated  with  type  3 
events.  It  is  clear  that  waves  near  the  lower  hybrid  fre¬ 
quency,  and  also  EMIC  waves  at  about  half  the  proton 
gyrofrequency,  can  be  generated  by  downgoing  keV  au¬ 


roral  electrons  [e.g.,  Chang  and  Coppi,  1981;  Vaivads  et 
al. ,  1998].  In  about  60  %  of  the  type  3  events  studied  by 
Norqvist  et  al.  [1998],  lower  hybrid  waves  seem  to  have 
higher  amplitude  than  all  other  types  of  waves  (type 
3(LH)  events).  In  about  20  %  of  the  events,  EMIC 
waves  are  the  most  intense  emissions  (type  3(EMIC)). 
For  the  remaining  20%  of  the  0+  heating  events,  broad¬ 
band  low-frequency  waves  have  the  highest  amplitudes. 
This  means  that  the  statistics  of  type  3  events  really 
gives  an  upper  limit  of  the  importance  of  EMIC  emis¬ 
sions  and  waves  near  the  lower  hybrid  frequency. 
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Figure  4.  Mean  0+  energy  as  a  function  of  energy  flux 
of  auroral  electrons  for  (a)  type  1  events  and  (b)  type  3 
events.  For  type  1  the  0+  energy  is  anticorrelated  with 
the  auroral  electrons  when  studied  on  timescales  of  a  few 
seconds,  while  for  type  3  the  0+  energy  is  correlated  with 
the  auroral  electrons.  (From  Andre  et  al.  [1998].) 


For  type  1  and  2  events,  broadband  low-frequency 
waves  dominate  the  ion  energization.  The  energy  trans¬ 
fer  from  the  waves  to  the  ions  can  be  modeled  and  it 
can  be  shown  that  the  wave  amplitudes  are  high  enough 
to  cause  the  observed  ion  energies  [e.g.,  Andre  et  al.y 
1998].  However,  the  origin  of  the  waves  is  not  under¬ 
stood.  One  possibility  is  that  the  braodband  waves 
contain  a  small  fraction  of  left-hand  polarized  electro¬ 
magnetic  component  (either  locally  generated  or  prop¬ 
agated  there),  which  can  gyroresonantly  energize  the 
ions  [Chang  et  al.7  1986]. 
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Alternatively,  essentially  electrostatic  waves  near  the 
ion  gyrofrequency  could  be  generated  by  a  field-aligned 
current  in  a  homogeneous  plasma  carried  by  a  drifting 
electron  distribution  [Kindel  and  Kennel ,  1971].  Us¬ 
ing  reasonable  assumptions,  Andre  et  a/.,  [1998]  showed 
that  the  observed  currents  are  not  high  enough  to  gener¬ 
ate  waves  in  this  way.  Also,  according  to  simple  theory, 
the  waves  that  would  be  generated  by  higher  current 
densities  might  not  be  as  broadband  as  the  observed 
emissions.  Field- aligned  electron  beams  with  higher 
energies  (up  to  hundreds  of  eV)  may  have  velocities 
matching  the  local  Alfven  speed,  When  these  beams 
have  high  enough  densities  they  may  generate  more 
electromagnetic  waves.  However,  as  for  the  electro¬ 
static  emissions,  these  waves  would  according  to  simple 
theory  not  be  as  broadband  as  the  observed  emissions. 

An  interesting  possibility  to  decrease  the  current 
needed  to  generate  waves,  and  to  obtain  more  broad¬ 
band  emissions,  is  to  introduce  gradients  perpendicular 
to  the  geomagnetic  field  in  the  model  [Ganguli,  1997; 
and  references  therein].  Laboratory  experiments  us¬ 
ing  parameters  relevant  for  space  plasmas  confirm  that 
gradients  in  perpendicular  electric  fields  decrease  the 
current  needed  for  wave  generation,  and  that  rather 
broadband  waves  can  be  obtained  [Amatucci  et  a/., 
1996;  Koepke  et  a/.,  1998;  and  references  contained 
therein].  A  possible  scenario  in  type  1  and  2  ion  ener¬ 
gization  regions  is  that  a  downward  propagating  Alfven 
wave  with  frequency  of  around  1  Hz  provides  a  gradi¬ 
ent  in  the  perpendicular  electric  field.  A  rather  small 
field-aligned  current  can  then  generate  fairly  broadband 
waves.  However,  this  scenario  is  hard  to  immediately 
confirm  from  observations.  It  is  not  clear  how  meaning¬ 
ful  it  is  to  split  the  measured  broadband  wave  spectrum 
into  one  part  much  below  the  ion  gyrofrequency  (the 
Alfven  wave  generated  somewhere  else,  providing  elec¬ 
tric  field  gradients)  and  a  more  broadband  part  around 
the  ion  gyrofrequency  (locally  generated  by  the  large 
scale  current,  in  the  presence  of  the  electric  field  gradi¬ 
ents).  The  latter  part  would  then  energize  the  ions.  A 
starting  point  for  a  detailed  investigation  could  be  the 
fact  that  the  ratio  of  the  electric  and  magnetic  fields 
of  the  broadband  waves  often  is  more  consistent  with 
Alfven  waves  at  the  lowest  frequencies,  as  compared 
with  the  more  electrostatic  emissions  at  higher  frequen¬ 
cies. 

6.  Summary 

More  than  one  ion  energization  mechanism  is  op¬ 
erating  in  the  terrestrial  auroral  region.  Using  Freja 
observations  from  altitudes  near  1700  km  in  the  auro¬ 
ral  region,  we  conclude  that  resonant  energization  by 
broadband  low-frequency  waves  is  the  most  important 
mechanism  for  the  energization  and  outflow  of  0+  ions. 


This  mechanism  dominates  the  0+  outflow  at  all  local 
times,  and  at  all  levels  of  magnetic  activity.  Resonant 
energization  by  waves  near  the  lower  hybrid  frequency, 
and  by  EMIC  waves,  occurs  in  the  pre-midnight  and 
midnight  sectors.  These  mechanisms  may  dominate 
specific  events,  but  give  lower  total  outflow  as  compared 
to  energization  by  broadband  low-frequency  waves. 
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Abstract.  It  is  shown  that  self-organized  criticality 
and  multiscale  intermittent  turbulence  characterize  the 
underlying  physics  of  the  substorm  behavior  of  the 
Earth’s  magnetotail.  The  model  incorporates  a  mul¬ 
titude  of  cross-scale  interactions  of  coherent  structures 
and  plasma  fluctuations.  Such  a  description  provides  a 
convenient  explanation  of  the  localized  and  sporadic  na¬ 
ture  of  the  reconnection  signatures  (bursty  bulk  flows) 
that  are  commonly  observed  in  the  tail  region  of  the 
magnetosphere  [Angelopoulos  et  al,  1996;  Lui,  1998, 
and  references  contained  therein].  The  onset  of  sub¬ 
storm  is  then  attributed  to  the  emergence  of  a  global 
nonclassical  nonlinear  instability  characterized  by  the 
enhanced  mixing  and  merging  of  the  coherent  struc¬ 
tures,  leading  to  the  phenomenon  of  self-organized  crit¬ 
icality  (SOC)  with  low-dimensional  chaotic  signatures 
and  fractal  behavior.  These  concepts  provide  a  new 
paradigm  for  the  understanding  of  the  ever-elusive  phe¬ 
nomenon  of  magnetic  substorms. 


1.  Synopsis 

In-situ  experimental  observations  [Angelopoulos 
et  al,  1996;  Lui,  1998,  and  references  contained  therein] 
and  simulations  [Drake,  1997;  Buchner ,  1998]  indicate 
that  the  magnetotail  dynamics  associated  with  sub¬ 
storms  are  intermittent  and  turbulent.  A  model  in¬ 
volving  the  generation,  dispersing,  and  merging  of  mul¬ 
tiscale  coherent  plasma  structures  and  fluctuations  is 
suggested  to  address  the  implications  of  such  observa¬ 
tions  [Chang,  1998]. 

The  merging  of  localized  coherent  structures  during 
such  a  dynamic  process  are  shown  to  be  responsible 
for  the  commonly  detected  sporadic  reconnection  sig¬ 
natures  and  associated  plasmoids  in  the  magnetotail 
and  contribute  to  the  origin  of  the  observed  “bursty 
bulk  flows” . 

The  commencement  of  a  substorm  is  relegated  to  the 
onset  of  a  global  nonclassical  nonlinear  instability,  sig¬ 
naled  by  the  enhanced  mixing  and  merging  of  coherent 
structures  and  other  plasma  fluctuations.  The  ensu¬ 
ing  substorm  dynamics  is  then  postulated  to  be  charac¬ 
terized  by  the  phenomenon  of  self-organized  criticality 
(SOC). 


Using  its  invariance  properties,  it  has  been  shown 
by  Chang  [1992],  that  systems  near  self-organized  crit¬ 
icality  will  generally  acquire  the  characteristics  of  low 
dimensionality  and  fractal  structures;  thus  providing  a 
rationale  for  the  conjectured  chaotic  behavior  of  sub¬ 
storm  dynamics  based  on  finite-dimensional  nonlinear 
dynamics  arguments  [Baker  et  al,  1990;  Vassiliadis  et 
al. ,  1990;  Klimas  et  al. ,  1992;  Sharma  et  ai,  1993; 
Klimas  et  ah,  1998]. 

This  leads  us  to  consider  the  fluctuation  spectra  for 
the  inhomogeneous,  anisotropic  intermittent  turbulent 
plasma  medium  of  the  magnetotail  during  substorms. 
The  concept  of  multifractals  is  then  introduced  and  re¬ 
lated  to  the  scaling  regions  of  the  spectra  for  coherent, 
MHD  and  kinetic  states.  Recent  observations  seem  to 
validate  the  above  hypothesis  [ffosino  et  ai,  1997;  Milo - 
vanov  et  ai ,  1996;  Zelenyi  et  ai,  1998]. 

We  will  conclude  by  making  a  few  suggestions  with 
regard  to  the  future  directions  of  observational,  lab¬ 
oratory  and  theoretical  research,  that  might  help  to 
elucidate  the  detailed  dynamics  of  this  proposed  new 
paradigm. 

2.  Plasma  Resonances  and  Coherent  Structures 

Measurements  of  reconnection  signatures  in  the 
“neutral  sheet”  region  of  the  Earth’s  magnetotail  in¬ 
dicate  that  most  of  the  individual  localized  merging 
processes  occur  at  intermediate  or  microscales.  On 
the  other  hand,  the  dimensions  of  the  full  dynamic 
domain  that  is  responsible  for  the  transferring  of  en¬ 
ergy  and  momentum  from  the  solar  wind  to  the  Earth’s 
magnetotail  generally  involve  time  and  spatial  scales 
much  larger  than  those  characterized  by  the  micro¬ 
scopic  plasma  parameters  such  as  the  ion  gyroradius, 
skin  depth,  Debye  length,  and  ion  cyclotron,  lower  hy¬ 
brid  or  plasma  frequencies.  The  transport  processes 
at  the  two  ends  of  the  dynamic  spectrum  can  involve 
characteristic  parameters  differ  by  orders  of  magnitude, 
suggesting  that  the  underlying  dynamics  of  the  magne¬ 
totail  is  intrinsically  multiscale. 

Thus,  we  adopt  a  middle  of  the  road  approach  in 
viewing  the  magnetotail  by  assuming  that  the  dynam¬ 
ics  of  the  plasma  medium  is  primarily  characterized  by 
the  basic  MHD  variables  with  an  anisotropic  pressure 
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tensor.  To  bring  in  some  of  the  possible  kinetic  effects, 
we  can,  e.g.,  relate  the  pressure  tensor  to  the  particle 
distribution  functions  /;(x,  v,2)  in  terms  of  the  appro¬ 
priate  moments. 

P  =  Em>  /(v-V)(v-V)/;(v)dv  (1) 

t  J 

where  /j(x,  v,£)  satisfy  the  Vlasov  equations. 

Standard  arguments  lead  to  the  following  coupled 
equations  of  induction  and  motion: 

dB/dt  =  (B  *  V)V  +  ...,  and  (2) 

p dV /dt—  (B  •  V)B  -f  . . . ,  (3) 

where  all  notations  are  customary  and  the  ellipsis  rep¬ 
resent  compressibility,  and  anisotropic  pressure  effects. 
Generally,  of  course,  dissipative  terms  must  also  be  in¬ 
cluded.  It  is  clear  from  above  that  one  of  the  wave 
modes  allowed  by  these  equations  is  the  shear  Alfven 
wave.  For  such  modes  to  propagate,  the  propagation 
vector  k  must  contain  a  field-aligned  component,  i.e., 
B  *  V  — y  ik  •  B  /  0.  However,  at  sites  where  the 
parallel  component  of  the  propagation  vector  vanishes, 
Arj|  =  k  *  B  =  0,  energies  are  localized  and  the  field 
lines  may  be  distorted  effortlessly.  These  singularities 
(points,  curves  or  surfaces)  at  which  =  0  are  called 
“shear  Alfven  resonances” .  As  it  will  be  demonstrated 
below,  the  existence  of  shear  Alfven  resonances  will  lead 
to  the  formation  of  nearly-nonpropagating  and  essen¬ 
tially  closed  macroscopic  magnetic  structures. 

We  now  consider  the  magnetic  field  structures  near 
the  Alfven  resonances.  Neglecting  the  pressure  effects, 
it  is  clear  from  Eqs.  (2),  (3)  that  the  forces  arise  from 
the  fluctuations  just  away  from  these  singular  (reso¬ 
nance)  sites,  i.e.,  JB  •  V,  will  tend  to  restore  the  field 
lines  towards  the  resonance  sites,  thereby  forming  es¬ 
sentially  closed  coherent  magnetic  structures.  In  the 
following  we  shall  consider  the  general  topology  of  such 
coherent  structures. 

For  an  ideal  MHD  system,  any  physically  acceptable 
magnetic  field  must  satisfy  V  •  B  =  0.  Also,  any  varia¬ 
tion  of  the  field  away  from  the  initial  value  must  satisfy 
the  constraint: 

dB/dt  =  Vx(vxB).  (4) 

Taylor  [1974]  demonstrated  that  Eq.  (4)  may  be  re¬ 
placed  by  an  infinite  set  of  integral  constraints  involving 
the  helicity  K ,  such  that 

I<  =  f  A  BdV  (5) 

Jv 

is  an  invariant  for  any  volume  V  enclosed  by  a  flux 
surface,  where  A  is  the  vector  potential.  It  can  be 
shown  that  as  the  system  relaxes  to  its  minimum  energy 
state  satisfying  the  helicity  conservation  constraints, 
the  magnetic  structure  will  be  in  a  force-free  state,  i.e., 

j  x  B  =  0.  (6) 


Figure  1.  Cross-sectional  view  of  sporadically  distributed 
flux  tubes  in  the  plasma  sheet. 


Taylor's  Conjecture  and  Coarse-Grained  Helicity .  Let 
us  now  consider  our  present  situation  at  hand.  We  are 
interested  in  the  more  realistic  situation  that  character¬ 
izes  the  dynamics  of  the  magnetotail  where  the  plasma 
is  slightly  dissipative  and  in  addition,  there  are  stochas¬ 
tic  macroscopic  (as  well  as  microscopic)  fluctuations. 
The  dissipation  and  magnetic  stochasticity  will  allow 
the  field  lines  to  merge,  mix,  and  break.  It  is  obvious 
that  it  no  longer  makes  sense  to  discuss  the  topology  of 
individual  field  lines.  Nevertheless,  it  was  suggested  by 
Taylor  [1974]  that  when  the  volume  integral  for  Eq.  (5) 
is  taken  over  the  “stochastic  region”,  the  coarse-grain 
averaged  helicity  will  be  essentially  conserved.  This  in¬ 
dicates  that  when  considering  the  stochastic  domain, 
the  average  magnetic  structure  in  a  relaxed  state  will 
again  be  essentially  force  free,  with  jxB  =  0,  where  j 
and  B  are  now  to  be  interpreted  as  the  mean  current 
and  magnetic  field,  respectively.  This  result  can  also  be 
arrived  at  using  the  clump  theory  of  MHD  turbulence 
[ Tetreault ,  1992,  and  references  contained  therein]. 

We  are,  of  course,  interested  in  the  magnetotail  at 
dynamic  states  that  are  far  from  equilibrium.  Thus, 
in  visualizing  the  relaxed  states  from  the  point  of  view 
of  the  Taylor’s  conjecture,  we  shall  consider  timescales 
such  that  “nearly  coherent”  magnetic  structures  are 
formed. 

Let  us  now  apply  these  concepts  to  the  sheared  mag¬ 
netic  field  geometries  that  are  generally  found  in  the 
“neutral  sheet”  region  of  the  magnetotail.  The  nearly 
force-free  condition  for  the  coarse-grain  averaged  co¬ 
herent  structures  would  then  orient  themselves  more- 
or-less  in  the  average  cross-tail  current  direction  in  the 
form  of  twisted  flux  tubes.  In  general,  there  will  be  a 
constellation  of  such  coherent  structures  immersed  in 
the  turbulent  plasma  medium,  Fig.  1. 
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Figure  2.  Cross-sectional  view  of  coherent  structures  at 
various  stages  of  merging,  (a)  Just  Prior  to  Merging,  (b) 
During  the  process  of  merging,  (c)  Relaxed  state  after  merg¬ 
ing. 


3.  Stochastic  Merging  of  Coherent  Structures 

As  the  coherent  structures  migrate  toward  each 
other,  they  will  merge  and  form  new  coherent  struc¬ 
tures.  Depending  on  the  polarities  and  intensities  of 
the  currents  that  orient  these  flux  tubes,  the  resulting 
coherent  structures  will  be  either  larger  or  smaller  than 
the  original  individual  structures.  The  final  states  of 
the  new  coherent  structures  will  again  be  essentially 
force- free  in  the  coarse-grained  sense.  As  these  new 
structures  are  generated,  new  MHD  fluctuations  are 
produced;  and  thereby  spontaneously  set  up  new  reso¬ 
nance  sites.  Thus,  an  interesting  scenario  of  intermit¬ 
tent  turbulent  mixing,  diffusing,  and  merging  sets  in. 

Let  us  consider  the  most  probable  situation  of 
merging,  i.e.,  the  merging  of  two  coherent  structures. 
Viewed  in  a  section  normal  to  the  average  direction  of 
the  cross- tail  current,  the  topologies  of  the  field  lines 
during  such  a  merging  process  mimic  that  is  generally 
considered  for  a  classical  magnetic  “reconnection”  pro¬ 
cess  [Fig.  2].  However,  we  note  that  this  localized  merg¬ 
ing  process  can  take  place  without  the  requirement  of 
B  =  0  and/or  the  existence  of  a  true  “neutral  line.”  In 
fact,  as  seen  above,  the  prerequisite  for  the  existence  of 
many  such  coherent  structures  as  well  as  the  sporadic 
merging  of  these  structures  is  the  existence  of  many 
“Alfven  resonance”  sites  with  =0.  This  occurs  when 
the  background  magnetic  field  is  three-dimensional  and 
nonzero  and  when  there  are  three-dimensional  macro¬ 
scopic  MHD  fluctuations. 

Thus,  we  suggest  that  as  a  spacecraft  flies  through 
the  neutral  sheet  region  of  the  magnetotail,  there  is  a  fi¬ 
nite  probability  for  the  instruments  on  the  spacecraft  to 
detect  classical-like  reconnection  signatures.  Such  sig¬ 


natures  can  be  detected  nearly  anywhere  in  the  plasma 
sheet,  but  more  probably  in  the  "neutral  sheet”  re¬ 
gion,  particularly  during  substorm  times.  The  dura¬ 
tion  of  interaction  of  these  observed  localized  merging 
processes  should  be  the  approximate  time  required  for 
the  new  relaxed  coherent  structures  to  emerge  and  in 
general,  would  be  rather  sporadic.  We  suggest  that 
these  are  the  origins  of  the  observed  “bursty  bulk  flows” 
[Angelopoulos  et  al. ,  1996;  Lui ,  1998;  Kivelson  and 
Kepko ,  private  communication]. 

Most  of  the  observed  localized  reconnection  signa¬ 
tures  to  date  seem  to  indicate  that  these  localized  merg¬ 
ing  processes  take  place  in  domain  sizes  comparable  to 
that  of  the  ion  gyroradius,  especially  during  substorm 
times.  Thus,  very  probably  most  of  these  processes 
will  be  influenced  by  microscopic  kinetic  effects.  Dur¬ 
ing  these  dynamic  processes,  the  ions  can  probably  be 
assumed  to  be  unmagnetized  and  the  electrons  fully 
magnetized  and  the  plasma  nearly  collisionless.  This, 
of  course,  would  lead  to  electron-induced  Hall  currents. 
Depending  on  the  underlying  magnetic  geometry  (since 
these  processes  can  occur  at  any  arbitrary  underlying 
magnetic  field  configuration),  the  relevant  kinetic  in¬ 
stability  that  can  initiate  the  localized  merging  (or  re¬ 
connection)  can  be  any  of  the  many  recently  suggested 
microscopic  instabilities  such  as  the  collisionless  tearing 
instability,  cross-field  two-stream  instability  [Lui,  1996], 
etc.  It  is  very  probable  that  the  nonlinear  state  of  merg¬ 
ing  for  each  of  these  localized  reconnections  again  en¬ 
tails  the  phenomenon  of  overlapping  resonances  [Galeev 
et  a/.,  1986].  (Now  these  resonances  will  arise  from 
the  localization  of  microscopic  fluctuations,  e.g.,  the 
whistler  resonances,  and  multiple  tearing  modes.) 

4.  The  Role  of  Self- Organized  Criticality  in  Sub¬ 
storm  Dynamics 

Under  favorable  conditions  (e.g.,  an  enhancement  of 
the  cross-tail  current  due  to  the  change  of  certain  global 
controlling  parameters  for  the  magnetotail),  the  state 
of  the  turbulence  discussed  in  the  previous  sections  may 
grow  (in  terms  of  more  plasma  fluctuations  and  larger 
coherent  structures). 

This  type  of  instability,  by  definition,  is  genuinely 
“nonlinear”.  For  the  onset  and  growth  of  a  classical 
nonlinear  instability,  there  generally  exists  a  prescribed 
minimum  finite  amplitude  of  disturbance  (measured, 
for  example,  by  the  root-mean-square  of  fluctuations) 
beyond  which  the  fluctuations  and  coherent  structures 
can  grow  [Fig.  3].  During  the  onset  of  a  substorm,  on 
the  other  hand,  it  has  been  recognized  [Chang,  1992] 
that  the  non-equilibrium  plasma  state  in  the  magneto¬ 
tail  is  near  criticality  (similar  to  the  critical  point  for 
equilibrium  liquid/gas  phase  transitions.)  At  such  a 
dynamic  state  (commonly  referred  to  as  a  state  of  self- 
organized  criticality,  SOC),  the  effect  of  the  fluctuations 
themselves  becomes  an  important  factor  in  determining 
the  critical  threshold  of  onset  and  the  instability  is  now 
both  “nonclassical”  and  "nonlinear”  [Fig.  3].  Such  an 
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Figure  3.  Marble  rolling  over  a  hill.  Schematic  depiction  of 
classical  and  nonclassical,  nonlinear  instabilities,  x*:  Clas¬ 
sical  threshold  for  nonlinear  instability. 


instability  is  generally  controlled  by  a  set  of  relevant 
global  parameters 

The  dynamics  of  a  system  at  self-organized  critical¬ 
ity  is  notoriously  difficult  to  handle  because  the  cor¬ 
relations  of  the  fluctuations  are  long-ranged  and  there 
exist  many  correlation  scales.  On  the  other  hand,  since 
the  correlations  are  extremely  long- ranged,  it  is  reason¬ 
able  to  expect  that  the  system  will  exhibit  some  sort 
of  invariance  under  scale  transformations  [Chang  et  a/., 
1992;  Chang ,  1992;  Chang ,  1998]. 

One  of  the  most  surprising  results  of  such  the¬ 
oretical  arguments  (based  on  the  ideas  of  the  dy¬ 
namic  renormalization-group)  is  that,  even  though  the 
magnetotail  system  is  infinite-dimensional  and  entails 
many,  many  physical  parameters,  its  dynamics  near 
self-organized  criticality  (i.e.,  at  the  onset  and  during 
magnetic  substorms)  can  actually  be  characterized  by 
a  small  number  of  relevant  physical  parameters  (low- 
dimensionality).  Thus,  the  system  can  be  approxi¬ 
mately  described  by  a  small  number  of  nonlinear  equa¬ 
tions  that  describe  the  evolution  of  these  relevant  phys¬ 
ical  parameters.  In  addition,  near  criticality,  these  low¬ 
dimensional  nonlinear  equations  will  generally  exhibit 
chaotic  (fractal)  behavior  [Chang,  1992]. 

In  several  recent  interesting  papers,  [Baker  et  a/., 
1990;  Klimas  et  a/.,  1991,  1992,  1998;  Baker ,  1998],  it 
has  been  suggested  that  certain  substorm  characteris¬ 
tics  in  terms  of  the  AL  time  series  could  be  modeled  by 
deterministic  chaos  of  simple  low-dimensional  dynamic 
equations.  These  results  suggest  that  the  global  mag- 
netospheric  system  may  be  “very  close”  to  forced  or 
self-organized  criticality. 

Vassiliadis  et  al.  [1990],  Shan  et  al.  [1991],  and 
Roberts  et  al  [1991]  applied  the  technique  of  Grass- 


Figure  4.  Fluctuation  spectrum  near  the  “neutral  sheet.” 

berger  and  Procaccia  [1983]  for  finite-dimensional  sys¬ 
tems  to  the  AE/AL  time  series  and  estimated  the  di¬ 
mensionally  of  the  magnetospheric  system  during  sub¬ 
storms.  Their  results  seemed  to  indicate  “fractals”  and 
low-dimensional  chaos. 

Separately,  Sharma  et  al.  [1993]  suggested  that 
the  dimensionality  of  the  magnetospheric  system  might 
be  obtained  using  the  singular  spectrum  analysis  con¬ 
structed  from  the  AE/AL  time  series.  Their  results 
also  seemed  to  indicate  low-dimensionality  and  chaos. 
These  results  appear  to  bestow  credence  to  the  hypoth¬ 
esis  that  substorm  dynamics  is  indeed  characterized  by 
the  phenomenon  of  self-organized  criticality. 

5.  Turbulent  Spectra  and  Multifractal  Behavior 

In  previous  sections,  we  demonstrated  that  the  dy¬ 
namics  of  the  magnetotail  is  characterized  by  multiscale 
intermittent  turbulence.  A  standard  technique  to  study 
the  behavior  of  such  type  of  turbulence  is  through  the 
properties  of  the  spectra  of  the  turbulent  fluctuations. 

For  example,  in  the  “neutral  sheet”  of  the  magneto¬ 
tail,  one  of  the  more  important  spectra  [Fig.  4]  to  con¬ 
sider  is  that  of  the  square  of  the  magnetic  fluctuations 
in  the  cross-tail  direction  (SB2).  We  expect  the  spec¬ 
tra  to  generally  exhibit  fractal  characteristics  (i.e.,  non- 
classical  slopes  with  discernible  deviations  from  those 
obtainable  by  naive  dimensional  arguments).  In  re¬ 
gions  where  the  fluctuations  and  merging  dimensions 
are  much  larger  than  that  of  the  local  ion  gyroradius, 
the  spectrum  is  expected  to  exhibit  two  distinguishable 
parts:  a  domain  characterized  by  the  larger  scale  coher¬ 
ent  structures  and  a  fractal  domain  characterized  by  the 
predominantly  MHD  fluctuations.  On  the  other  hand, 
in  regions  within  the  narrow  cross-tail  current  sheet, 
we  expect  the  spectra  to  exhibit  at  least  three  distin¬ 
guishable  parts:  a  domain  that  contains  predominantly 
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large  scale  coherent  structures,  an  MHD  fractal  domain 
and  a  kinetic  fractal  regime  whose  fractal  dimension (s) 
generally  depends  on  the  type(s)  of  microscopic  fluctu¬ 
ations  and  microinstabilities  that  are  relevant  for  the 
merging  and  diffusion  processes.  Such  type  of  fluctua¬ 
tion  spectra  has  been  recently  observed  [Hosino  et  a/., 
1994;  Milovanov  et  a/.,  1996;  Zelenyi  et  a/.,  1998;  Lui , 
1998;  and  references  contained  therein.]. 

The  shapes  (slopes)  of  these  spectra  in  the  distant- 
tail  region  have  been  compared  with  results  based  on 
theoretical  scaling  ideas  involving  fractal  dimensions 
[Milovanov  et  al. ,  1996;  Zelenyi  et  al ,  1998].  The  dif¬ 
ference  of  slopes  of  the  various  domains  of  an  individ¬ 
ual  spectrum  indicates  that  the  scaling  (fractal)  behav¬ 
ior  of  each  domain  belongs  to  a  different  “universality 
class”.  Such  type  of  change  of  scaling  behavior  from 
one  universality  class  to  another  is  called  “symmetry 
breaking” .  In  addition  to  the  scaling  properties  of  indi¬ 
vidual  discernible  domains,  there  are  also  intermediate 
regimes  whose  fractal  properties  are  much  more  compli¬ 
cated  (as  indicated  by  the  circled  region  of  Fig.  4).  The 
scaling  laws  for  these  regions  are  generally  expected  to 
exhibit  multiple- power  or  other  nonlinear  characteris¬ 
tics  [Chang  et  al .,  1992]. 


6.  Recapitulation  and  Suggested  New  Research 
Directions  of  Magnetotail  Dynamics  during  the 
Next  Millennium 

In  summary,  we  have  introduced  a  multiscale  in¬ 
termittent  turbulence  model  for  the  dynamics  of  the 
magnetotail  [Chang,  1998].  The  theory  is  based  on 
the  overlapping  resonances  of  plasma  fluctuations.  It 
provides  a  physical  picture  of  sporadic  and  localized 
merging  of  coherent  magnetic  structures  of  varied  sizes. 
Such  a  picture  seems  to  depict  the  observational  prop¬ 
erties  of  “bursty  bulk  flows”  (sporadic  localized  recon¬ 
nections)  in  the  magnetotail  [Angelopoulos  et  al .,  1996; 
Lui,  1998].  In  this  picture,  the  onset  of  substorm  is 
due  to  a  global  nonclassical  nonlinear  instability  and 
the  dynamics  of  the  magnetotail  during  the  evolution 
of  the  substorm  is  characterized  by  the  phenomenon  of 
forced  or  self-organized  criticality  [Chang,  1998]. 

The  consequence  of  this  is  the  prediction  of  multi¬ 
fractal  characteristics  of  the  fluctuation  spectra  [Hosino 
et  al,  1994;  Milovanov  et  al.,  1996;  Zelenyi  et  al.,  1998; 
Chang,  1992;  Chang,  1998;  and  references  contained 
therein]  and  the  dynamics  of  the  magnetotail  behaves 
essentially  as  a  low  dimensional  system.  This  conclu¬ 
sion  seems  to  agree  with  the  results  of  some  of  the  recent 
nonlinear  dynamics  calculations  [Baker  et  al.,  1990;  Kli- 
mas  et  al.,  1991,  1992,  1998;  Baker,  1998]. 

If  we  accept  this  new  paradigm  of  magnetotail  dy¬ 
namics,  particularly  during  the  periods  of  substorm 
evolution,  it  immediately  indicates  several  directions  of 
observational,  numerical  and  laboratory  investigations 
for  the  21st  millennium. 


6.1.  Observations 

In  the  area  of  observations,  it  is  obvious  that  more 
information  in  the  direction  of  the  fluctuation  spectra 
will  be  needed.  With  the  limitation  of  single  spacecraft 
observations,  the  spectra  that  can  be  obtained  based 
on  present  day  observational  capabilities  will  necessar¬ 
ily  be  restricted  within  the  time  domain.  To  convert 
these  to  the  k  spectra,  some  assumptions  will  have  to  be 
made.  In  the  MHD  domain,  this  is  somewhat  straight¬ 
forward  since  the  wave  speed  is  essentially  known.  For 
the  kinetic  region,  however,  some  information  with  re¬ 
gard  to  the  wave  modes  must  be  available.  In  the  fu¬ 
ture,  when  sophisticated  coordinated  multi-spacecraft 
observations  are  available,  it  might  be  useful  to  mea¬ 
sure  the  multipoint  correlation  spectra  as  well. 

It  will  also  be  useful  if  more  statistically  averaged 
information  is  obtained  through  observations.  Since 
the  turbulence  is  anisotropic  and  inhomogeneous,  direc¬ 
tional  as  well  as  spatial  variations  are  needed  as  inputs 
in  the  development  of  specific  theoretical  calculations. 
After  all,  we  still  need  to  determine  what  is  (are)  ex¬ 
actly  the  global  instability  (instabilities)  that  is  (are) 
responsible  for  the  onset  of  substorms.  There  may  be 
several  different  types  of  instabilities  that  can  lead  to 
the  onset  of  substorms.  There  may  also  be  several  dif¬ 
ferent  classes  of  substorms,  depending  on  the  input  of 
the  global  parameters;  thereby  leading  to  different  char¬ 
acteristics  of  the  fluctuation  spectra. 

6.2.  Numerical  simulations 

Recent  advances  in  global  simulations  have  given 
some  insights  to  the  nature  and  characteristics  of  the 
global  instability  (instabilities)  that  are  associated  with 
the  onset  and  evolution  of  substorms.  However,  since 
the  global  simulations  generally  depend  on  ad  hoc  as¬ 
sumptions  of  anomalous  resistivity  or  numerically  gen¬ 
erated  dissipation,  these  numerical  results  cannot  pro¬ 
vide  any  hint  to  the  nature  of  the  true  connection  be¬ 
tween  the  localized  (and  at  least  partially  kinetic)  dissi¬ 
pative  processes  and  the  apparent  dissipation  that  are 
required  in  most  of  the  global  simulation  calculations. 
Recent  numerical  simulations  of  reconnection  processes 
have  made  numerous  advances  in  the  understanding  of 
localized  reconnections.  These  results  are  particularly 
useful  in  the  interpretation  of  the  kinetic  portion  of  the 
fluctuation  spectra.  However,  such  information  alone 
will  not  be  sufficient  to  bridge  the  gap  of  the  global  pro¬ 
cesses  and  the  localized  dissipations.  An  entirely  differ¬ 
ent  class  of  numerical  simulations  must  be  performed. 
The  simulations  must  address  the  “large”  magnetic 
Reynolds  number  problem,  and  the  resulting  turbulence 
must  be  anisotropic,  inhomogeneous,  and  intermittent 
such  that  they  are  capable  of  resolving  the  generation, 
mixing,  and  merging  of  coherent  magnetic  structures 
as  well  as  the  interaction  of  these  effects  with  macro¬ 
scopic  and  microscopic  fluctuations  (including  particle 
acceleration  and  diffusion)  during  substorms. 


6.3.  Laboratory  simulations 

There  have  been  interesting  laboratory  simulations 
on  magnetic  reconnections.  For  example,  there  is  the 
admirable  work  of  Yamada  et  al  [1990]  on  the  merging 
of  MHD  flux  structures,  and  the  equally  fine  observa¬ 
tions  on  the  whistler  modes  generated  during  electron 
MHD  reconnection  processes  by  Stenzel  et  al  [1996]. 
Gekelman  and  his  collaborators  [1997]  have  been  mak¬ 
ing  very  interesting  and  accurate  studies  of  the  dynam¬ 
ics  and  merging  of  kinetic  Alfven  resonance  cones.  It 
will  be  instructive  (particularly  to  the  theorists)  when 
the  work  of  Gekelman  et  al  and  Stenzel  are  extended 
to  the  nonlinear  regime  of  overlapping  kinetic  Alfven 
and  whistler  resonances.  Ultimately  it  will  be  useful 
in  developing  larger  plasma  chambers  so  that  bonafide 
MHD  intermittent  turbulence  may  be  conducted. 

7.  Addendum 

In  an  interesting  recent  paper,  Chapman  et  al  [1998] 
presented  the  simulation  results  of  a  simple  avalanche 
model,  that  displayed  many  of  the  characteristic  fea¬ 
tures  (such  as  self-organized  criticality  and  global  in¬ 
stability,  etc.)  of  magnetospheric  activity  that  was  ad¬ 
vocated  in  this  treatise. 
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Abstract.  Recent  in  situ  observations  have  revealed 
novel  features  in  the  polar  wind.  Measurements  be¬ 
tween  5000  and  9000  km  altitude  by  the  Akebono  satel¬ 
lite  indicate  that  both  H+  and  0+  ions  have  remarkably 
higher  outflow  velocities  in  the  sunlit  region  than  on  the 
nightside.  Electrons  also  display  an  asymmetric  behav¬ 
ior:  the  dayside  difference  in  energy  spread,  greater 
for  upward-moving  than  downward-moving  electrons, 
is  absent  on  the  nightside.  Our  previous  calcula¬ 
tions  based  on  a  self-consistent  hybrid  model  [Tam 
et  a/.,  1995b]  have  demonstrated  the  significance  of  the 
anisotropic  kinetic  effects  of  photoelectrons  in  the  day- 
side  polar  outflow.  Results  generated  from  the  model 
agree  with  various  qualitative  features  observed  in  the 
sunlit  polar  region  by  the  Akebono  satellite.  Our  re¬ 
cent  work  has  extended  the  application  of  the  model  to 
the  night-side  polar  outflow,  in  which  photoelectrons 
are  practically  absent.  By  comparing  the  daytime  and 
night-time  results  of  our  model,  we  demonstrate  the 
anisotropic  kinetic  photoelectron  effects  on  the  polar 
outflow.  In  particular,  the  presence/ absence  of  these 
suprathermal  electron  effects  is  responsible  for  the  polar 
wind  day-night  asymmetries  observed  by  the  Akebono 
satellite. 

1.  Motivation 

The  existence  of  the  polar  wind,  an  outflow  of  plasma 
along  the  open  magnetic  field  lines  emanating  from  the 
polar  region  of  the  ionosphere,  was  first  proposed  by 
Axford  [1968]  and  Banks  and  Holzer  [1968].  These  early 
studies  recognized  the  ambipolar  electric  field  as  one  of 
the  mechanisms  governing  the  plasma  outflow.  Because 
the  polar  cap,  in  general,  is  a  relatively  quiescent  region, 
the  ambipolar  effect  is  a  major  contribution  to  the  elec¬ 
tric  field  in  the  “classical”  polar  wind,  the  steady-state, 
quasi- neutral,  current-free  outflow  of  plasma. 

The  ambipolar  field  itself  exists  self-consistently  with 
the  background  plasma.  It  is  therefore  influenced  by 
other  mechanisms  that  need  to  be  included  in  the  dy¬ 
namics  of  the  particles.  For  example,  the  geomagnetic 
field,  which  decreases  with  altitude,  gives  rise  to  the 
mirror  force  that  changes  the  particles’  pitch  angles. 
Coulomb  interactions  among  all  the  species  lead  to  en¬ 
ergy  exchange  and  pitch  angle  diffusion.  These  effects 
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are  essential  to  the  dynamics  of  the  particles,  and  there¬ 
fore,  can  affect  the  polar  wind  electric  field. 

Recent  observations  have  suggested  that  another  ef¬ 
fect,  photoelectron  populations  generated  in  the  sunlit 
ionosphere,  can  alter  the  polar  wind  significantly.  Our 
study  on  the  photoelectron-driven  polar  wind  is  mo¬ 
tivated  by  these  increasingly  convincing  experimental 
indications. 

Early  polar  cap  measurements  obtained  by  the 
ISIS-1  satellite  showed  evidence  of  “anomalous”  field- 
aligned  photoelectron  fluxes  in  both  upward  and  down¬ 
ward  directions,  where  the  downgoing  (return)  fluxes 
were  considerably  smaller  than  the  outgoing  fluxes 
above  a  certain  energy  [Winningham  and  Heikkilay 
1974].  Such  non-thermal  features  were  confirmed  by 
the  DE-1  and  -2  satellites  [Winningham  and  Gurgiolo , 
1982]:  outgoing  field- aligned  electron  fluxes  in  the  pho¬ 
toelectron  energy  range  were  observed  by  the  HAPI 
(High  Altitude  Plasma  Instrument)  on  DE-1  and  the 
LAPI  (Low  Altitude  Plasma  Instrument)  on  DE-2;  ev¬ 
idence  of  downstreaming  electron  fluxes  was  also  found 
in  the  low- altitude  distribution  measured  by  the  LAPI. 
These  fluxes  are  considered  anomalous  because  their  ex¬ 
istence  cannot  be  related  to  the  idea  of  thermal  conduc¬ 
tivity  and  temperature  gradient  in  classical  fluid  theo¬ 
ries.  Similar  to  the  ISIS-1  measurements,  the  return 
fluxes  observed  by  DE-2  were  comparable  to  the  outgo¬ 
ing  fluxes  below  some  truncation  energy,  but  consider¬ 
ably  smaller  above  that.  As  suggested  by  Winningham 
and  Gurgiolo  [1982],  the  existence  of  such  downstream¬ 
ing  fluxes  may  be  due  to  reflection  of  electrons  by  the 
ambipolar  electric  field  along  the  geomagnetic  field  line 
above  the  satellite.  The  truncation  energy,  obtained  by 
comparing  the  outgoing  and  the  return  electron  fluxes, 
would  thus  provide  an  estimate  for  the  potential  drop 
due  to  the  electric  field.  These  authors  observed  that 
this  truncation  energy  ranged  from  5  to  60  eV,  and 
thus  were  able  to  deduce  the  magnitude  of  the  poten¬ 
tial  drop  above  the  altitude  of  the  satellite  (~  500  km). 
Unfortunately,  existing  classical  polar  wind  theories 
can  only  account  for  a  much  smaller  potential  drop 
[ Ganguli ,  1996,  and  references  therein].  Winningham 
and  Gurgiolo  [1982]  also  pointed  out  that  variation  of 
the  truncation  energy  was  due  to  changes  in  the  solar 
zenith  angle  at  the  production  layer  below  the  satellite. 
The  solar  zenith  angle  is  related  to  the  photoionization 
rate,  which  itself  is  related  to  the  local  ionospheric  pho- 
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toelectron  density  [J asperse,  1981].  These  observations 
therefore  imply  a  relationship  between  the  local  pho¬ 
toelectron  density  below  the  satellite  and  the  potential 
drop  along  the  field  line  above  it,  and  are  consistent 
with  the  idea  that  the  photoelectrons  may  significantly 
affect  the  ambipolar  electric  field. 

While  the  observations  discussed  above  have  implied 
that  photoelectrons  may  contribute  to  the  dynamics  of 
the  polar  wind,  more  recent  evidence  indicates  that  the 
polar  wind  characteristics  themselves  are  affected  by 
the  photoelectrons.  In  situ  measurements  by  the  Ake- 
bono  satellite  have  revealed  novel  features  in  the  polar 
wind:  day-night  asymmetries  in  the  ion  and  electron 
features.  The  most  dramatic  are  asymmetries  in  the 
ion  outflow  velocities  [Abe  et  ai,  1993b]:  satellite  data 
between  5000  and  9000  km  altitude  have  indicated  re¬ 
markably  higher  outflow  velocities  for  the  major  ion 
species,  H+  and  0+,  in  the  sunlit  region  than  on  the 
nightside.  For  example,  the  H+  velocity  (u/,)  was  found 
to  be  about  12  km/s  on  the  dayside,  but  only  about 
5  km/s  on  the  nightside.  Similarly,  the  0+  velocity 
(u0)  in  the  sunlit  region  (~  7  km/s)  is  about  twice  that 
in  the  midnight  sector  (~  3  km/s).  A  day-night  asym¬ 
metry  was  also  observed  in  the  electron  behavior.  Elec¬ 
trons  were  distinguished  according  to  their  velocities 
along  the  geomagnetic  field  line.  On  the  dayside,  it  was 
found  that  the  temperature  of  the  upstreaming  popu¬ 
lation  is  greater  than  that  of  the  downstreaming  popu¬ 
lation,  i.e.,  Te>up  >  Te}down ,  indicative  of  an  upwardly 
directed  heat  flux  [Yau  et  ai ,  1995].  On  the  nightside, 
in  contrast,  no  such  up-down  anisotropy  was  observed 
[Abe  et  ai,  1996]. 

Besides  the  day-night  asymmetries,  Akebono  mea¬ 
surements  between  5000  and  9000  km  altitude  have 
also  revealed  other  sometimes  unexpected  ion  transport 
properties  in  the  polar  region  [Abe  et  ai,  1993b].  For 
example,  0+  was  most  often  found  to  be  dominant  over 
H+  as  the  major  ion  species,  contrary  to  the  traditional 
belief  that  very  few  0+  ions  are  able  to  overcome  the 
gravitational  force  and  escape  to  such  high  altitudes  due 
to  their  heavier  mass.  The  measured  outflow  velocities 
for  both  the  H+  and  0+  ions  in  general  increase  mono- 
tonically  with  altitude,  and  the  flows  for  both  species 
are  supersonic  at  high  altitudes.  In  fact,  the  measured 
0+  outflow  velocities  (see  above)  are  much  larger  than 
the  values  expected  by  classical  polar  wind  models  [e.g. 
Schunk  and  Watkins ,  1981,  1982;  Blelly  and  Schunk , 
1993].  All  these  ion  outflow  characteristics,  particularly 
the  enhanced  ion  outflow  velocities,  suggest  a  higher 
ambipolar  electric  field  than  that  predicted  by  classi¬ 
cal  polar  wind  models  [ Ganguli ,  1996,  and  references 
therein],  and  are  consistent  with  the  values  of  the  field- 
aligned  potential  drop  deduced  by  Winningham  and 
Gurgiolo  [1982]  based  on  the  DE-2  measurements. 

Because  of  the  marked  day-night  asymmetries  ob¬ 
served  in  several  characteristics  of  the  polar  wind,  and 
the  fact  that  photoelectrons  exist  primarily  in  the  sunlit 
ionosphere,  they  are  the  natural  candidate  to  account 
for  the  day-night  asymmetries.  Indeed,  collisionless  ki¬ 


netic  calculations  by  Lemaire  [1972]  showed  that  escap¬ 
ing  photoelectrons  may  enhance  the  electric  field  and 
increase  the  ion  outflow  velocities  in  the  polar  wind. 
Photoelectrons,  therefore,  may  provide  a  possible  ex¬ 
planation  for  the  observations  of  both  sets  of  satellites: 
the  magnitude  of  the  ambipolar  electric  field  deduced 
from  the  DE-2  measurements,  and  the  day-night  asym¬ 
metries  and  enhanced  ion  outflow  velocities  observed 
by  the  Akebono  satellite.  Because  Coulomb  collisions 
may  also  influence  the  dynamics  of  the  photoelectrons, 
for  example,  by  transferring  their  energy  to  other  par-  ^ 
tide  components  in  the  polar  wind,  and  thereby  re¬ 
ducing  the  escaping  photoelectron  flux,  collisional  ef¬ 
fects  should  also  be  taken  into  account  in  determining 
the  impact  of  photoelectrons  on  the  electric  field.  Our 
goal,  therefore,  is  to  address  these  observations  by  in¬ 
corporating  the  complete  photoelectron  physics  into  a  - 
self-consistent,  global  description  of  the  polar  wind. 

2.  Photoelectrons,  Energy  Fluxes  and  Electric 
Field 

The  global  kinetic  collisional  physics  of  suprather- 
mal  electrons  in  a  steady-state  space  plasma  outflow 
was  first  considered  by  Scudder  and  Olbert  [1979]  in 
their  study  of  the  solar  wind  halo  electrons.  These  au¬ 
thors  related  the  anomalous  field-aligned  electron  heat 
fluxes  observed  in  the  solar  wind  to  the  non-local  na¬ 
ture  of  the  electron  distributions,  and  demonstrated  the 
formation  of  such  non-thermal  features  using  a  simpli¬ 
fied  collisional  operator.  They  also  suggested  that  these 
suprathermal  electrons,  through  their  anomalous  con¬ 
tribution  to  the  energy  flux,  may  significantly  increase 
the  ambipolar  electric  field  along  the  magnetic  field 
lines,  thereby  “driving”  the  solar  wind  [Olbert,  1982]. 

An  analogous  situation  exists  for  the  dayside, 
photoelectron-driven  polar  wind.  It  has  been  shown  by 
Yasseen  et  ai  [1989]  that  the  polar  wind  photoelectrons 
can  give  rise  to  the  non-thermal  distributions  observed 
by  the  DE  satellites  [Winningham  and  Gurgiolo ,  1982]. 

The  effect  on  the  polar  wind  due  to  the  energy  fluxes  as¬ 
sociated  with  these  photoelectrons  has  been  examined 
by  Tam  et  ai  [1995a],  who  concluded  that  such  anoma¬ 
lous  electron  energy  fluxes  may  significantly  increase  ^ 

the  ambipolar  electric  field.  Because  photoelectrons  ex¬ 
ist  primarily  in  the  sunlit  ionosphere,  they  enhances  the 
dayside  ambipolar  electric  field,  thereby  increasing  the  ^ 

ion  outflow  velocities  on  the  dayside.  Photoelectrons, 
with  their  associated  energy  fluxes,  can  therefore  pro¬ 
vide  not  only  a  mechanism  for  the  enhanced  ion  outflow 
velocities  observed  by  the  Akebono  satellite  [Abe  et  ai, 

1993a,  1993b],  but  also  the  explanation  for  the  observed 
day-night  asymmetric  ion  and  electron  features  in  the 
polar  wind  [Abe  et  ai,  1993b,  1996;  Yau  et  ai,  1995]. 

Energetic  suprathermal  electrons  in  the  polar  wind 
or  in  other  ionospheric/magnetospheric  settings  have 
been  considered  by  various  authors.  For  example,  ki¬ 
netic  collisional  calculations  by  Khazanov  et  ai  [1993] 
have  examined  the  role  of  photoelectrons  on  plasma- 
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spheric  refilling.  Collisionless  kinetic  calculations  by 
Lemaire  [1972]  have  shown  that  escaping  photoelec¬ 
trons  may  increase  the  ion  outflow  velocities  in  the  po¬ 
lar  wind.  Collisionless  kinetic  calculations  by  Barakat 
and  Schunk  [1984]  and  generalized  semi-kinetic  (GSK) 
calculations  by  Ho  et  ah  [1992]  have  examined  the  im¬ 
pact  of  hot  magnet ospheric  electrons,  and  concluded 
that  such  particles  may  also  increase  the  ion  outflow 
velocities. 

We  should  add  that  other  mechanisms  besides 
suprathermal  electron  effects  may  also  be  proposed  as 
alternative  explanations  for  the  enhanced  ion  velocities. 
Parallel  ion  acceleration  driven  by  E  x  B  convection 
was  considered  by  Cladis  [1986],  and  shown  to  signifi¬ 
cantly  energize  0+  ions  escaping  to  the  polar  magneto¬ 
sphere.  This  force  can  also  be  seen  as  a  centrifugal  force 
in  the  convecting  frame  of  reference,  and  was  included 
in  this  form  in  the  time-dependent,  GSK  model  devel¬ 
oped  by  Horwitz  et  al  [1994].  These  mechanisms,  in¬ 
cluding  the  suprathermal  electron  effects,  have  recently 
been  reviewed  by  Ganguli  [1996]. 

Recently,  a  self-consistent  hybrid  model  has  been 
developed  by  Tam  et  al.  [1995b]  to  take  into  account 
the  global  kinetic  collisional  nature  of  the  polar  wind 
physics  introduced  by  the  photoelectrons.  The  details 
and  advantages  of  the  model  have  recently  been  dis¬ 
cussed  by  Tam  et  al.  [1998].  Our  results  in  this  pa¬ 
per  are  based  on  the  model.  The  model  represents 
two  breakthroughs  in  polar  wind  theoretical  model¬ 
ing.  One,  it  was  the  first  to  successfully  incorporate 
the  global  kinetic  collisional  photoelectron  effects  into 
a  self-consistent  polar  wind  description.  Two,  due  to  its 
kinetic  treatment  of  the  ions,  the  model  was  the  first  to 
generate  self-consistent  global  polar  wind  calculations 
whose  solutions  span  continuously  from  a  collisionless 
subsonic  regime  at  low  altitudes  to  a  collisional  super¬ 
sonic  regime  at  high  altitudes. 

The  model  is  hybrid  in  that  it  consists  of  a  kinetic 
and  a  fluid  component.  Photoelectrons  (treated  as  test 
particles  because  of  their  low  relative  density)  and  both 
the  H+  and  0+  ions  are  described  using  a  global  ki¬ 
netic  collisional  approach  while  thermal  electron  prop¬ 
erties  (density,  drift  velocity,  and  temperature)  are  de¬ 
termined  from  a  simpler,  fluid  approach  that  also  cal¬ 
culates  the  self-consistent  ambipolar  electric  field.  Be¬ 
cause  of  its  treatment  of  the  thermal  electrons,  the 
model  should  be  distinguished  from  traditional  hybrid 
approaches  where  electrons  are  treated  as  a  massless 
neutralizing  fluid.  The  model  is  based  on  an  iterative 
scheme  combining  the  kinetic  and  fluid  calculations, 
that  should  converge  to  physically  meaningful  solutions. 

Let  us  specify  some  criteria  that  will  enable  us  to  de¬ 
fine  more  precisely  our  study  of  the  polar  wind.  First, 
we  will  consider  the  polar  wind  only  at  altitudes  above 
500  km  (which  corresponds  roughly  to  the  polar  or¬ 
bits  of  DE-2).  At  such  altitudes,  neutral  densities  are 
low  enough  to  neglect  the  “chemical”  reactions  such  as 
photoionization,  recombination,  etc.  Second,  the  mag¬ 
nitude  of  the  geomagnetic  field  is  such  that  the  gyration 


period  and  Larmor  radius,  for  all  particle  species,  are 
much  smaller  than  any  relevant  time  or  length  scales. 
We  can  therefore  use  the  guiding  center  approximation. 
Third,  the  gradients  of  the  geomagnetic  field  are  such 
that  only  transport  along  the  geomagnetic  field  line  is 
important.  The  time-dependent  distribution  function 
f(t}s,v ||,u_l)  for  a  given  particle  species  is  therefore 
governed  by  the  following  collisional  gyrokinetic  equa¬ 
tion: 
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where  s  is  the  distance  along  the  magnetic  field  line  B ,  q 
and  m  are  the  algebraic  electric  charge  and  mass  of  the 
species  respectively,  E\\  is  the  field- aligned  electric  field, 
g  is  the  gravitational  acceleration,  B'  =  dB/ds ,  8f/St 
represents  the  rate  of  change  of  the  distribution  func¬ 
tion  due  to  collisions,  and  C  is  a  collisional  operator  for 
Coulomb  interaction,  which  is  the  dominant  type  of  col¬ 
lision  above  300  km  altitude.  Note  that  the  operator  C 
for  a  given  particle  species  describes  not  only  Coulomb 
collisions  with  other  species,  but  also  those  among  the 
same  species  itself.  Equation  (1)  thus  includes  all  the 
major  forces  a  particle  experiences  as  it  travels  along 
the  geomagnetic  field  line:  gravitational  force,  field- 
aligned  electric  force,  mirror  force,  and  forces  that  are 
due  to  Coulomb  collisions,  including  those  with  the 
same  species.  The  time-independent  version  of  Eq.  (1) 
is  the  governing  equation  in  the  kinetic  component  of 
our  model. 

3.  Day-Night  Asymmetries  in  the  Polar  Outflow 

The  main  goal  of  this  paper  is  to  demonstrate  the  ob¬ 
served  day-night  asymmetry  of  the  polar  outflow  based 
on  our  study  of  the  anisotropic  kinetic  photoelectron  ef¬ 
fects.  Application  of  the  self-consistent  hybrid  model  to 
the  sunlit  polar  region  has  generated  results  [Tam  et  a/., 
1998]  that  are  qualitatively  and  quantitatively  consis¬ 
tent  with  polar  wind  observations  [Abe  et  a/.,  1993b; 
Yau  et  al. ,  1995].  In  this  study,  we  characterize  the 
night-time  polar  wind  conditions  with  the  absence  of 
photoelectrons.  Thus,  we  generate  a  night-time  polar 
outflow  solution  by  using  the  same  boundary  conditions 
as  those  for  the  daytime,  except  without  the  presence 
of  photoelectrons,  i.e.  ns  =  0. 

The  density  profiles  of  all  the  species  in  the  day¬ 
time  polar  wind  solution  is  shown  in  Fig.  1.  Note  that 
the  photoelectron  density  ns  is  small  compared  with 
the  thermal  electron  density  ne.  Our  test-particle  ap¬ 
proach  for  photoelectrons  is  therefore  justified.  The 
small  photoelectron  density,  however,  does  not  imply 
that  these  suprathermal  particles  are  insignificant  to 
the  overall  polar  outflow  dynamics.  Because  of  the  re¬ 
lation  between  the  electron  energy  flux  (or  heat  flux) 
and  the  ambipolar  electric  field  [Tam  et  a/.,  1995a],  one 
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Figure  1.  Density  profiles  for  0+  (n0),  H+  («>,),  thermal 
electron  (ne),  and  photoelectron  (n3)  in  the  dayside  polar 
wind  solution.  Note  that  ne  and  n0  are  almost  equal  because 
of  quasi-neutrality.  Their  lines  virtually  overlap  in  the  plot. 


Figure  2.  A  typical  normalized  distributions  for  the  com¬ 
bined  electron  population  (thermal  and  photoelectrons)  in 
the  dayside  polar  wind  (solid  line).  Dot-dashed  line  is  for 
the  thermal  electron  distribution;  dashed  line  is  for  the  pho- 
toelectrons.  The  distributions  were  calculated  at  7500  km 
altitude.  Te  is  the  temperature  of  the  thermal  electron  pop¬ 
ulation,  and  me  is  the  electron  mass. 


has  to  examine  the  photoelectron  heat  flux  contribution 
as  well  in  order  to  determine  whether  the  suprathermal 
population  can  significantly  affect  the  outflow  dynam¬ 
ics.  Fig.  2  shows  a  typical  normalized  distribution  for 
the  combined  electron  population  in  our  dayside  so¬ 
lution.  The  contributions  from  the  thermal  and  pho¬ 
toelectrons  are  also  plotted  in  the  figure.  Note  that 
the  tail  portion  of  the  distribution  is  dominated  by  the 
photoelectron  population.  Such  a  distribution  suggests 
that  the  majority  of  the  electron  heat  flux  is  carried  by 
the  photoelectrons.  By  comparing  the  heat  flux  contri¬ 
bution  due  to  the  thermal  electron  with  that  carried  by 
the  photoelectrons,  we  find  that  the  heat  fluxes  carried 
by  the  two  electron  components  are  in  opposite  direc¬ 
tions:  downward  for  the  thermal  but  upward  for  the 
suprathermal,  as  shown  in  the  top  panel  of  Fig.  3.  We 
also  find  that  the  photoelectron  heat  flux  is  much  larger 
in  magnitude  than  its  thermal  counterpart,  as  shown  in 


Figure  3.  Top  panel:  heat  fluxes  carried  by  the  thermal 
electrons  (Qe, thermal)  and  photoelectrons  (Qe, suprathermal) 
in  the  dayside  polar  wind  solution.  Note  that  the  heat  flux 
contribution  by  the  thermal  electrons  is  in  the  downward 
direction  (negative  sign)  while  that  by  the  photoelectrons  is 
upwardly  directed  (positive  sign).  Bottom  panel:  the  ratio 
of  the  magnitudes  of  the  two  heat  fluxes. 


the  bottom  panel  of  Fig.  3.  The  two  heat  flux  compo¬ 
nents  combine  to  give  the  total  electron  heat  flux  whose 
direction  is  dictated  by  the  photoelectron  contribution, 
i.e.  upward.  The  upwardly  directed  total  electron  heat 
flux  is  consistent  with  the  data  from  Yau  et  al  [1995] 
and  the  results  by  Tam  et  al  [1995a,  1995b], 

Due  to  the  upwardly  directed  heat  flux  associated 
with  the  kinetic  nature  of  photoelectrons,  we  expect 
the  self-consistent  ambipolar  electric  field  to  be  larger  in 
the  dayside  polar  wind.  That  would  give  rise  to  a  larger 
potential  difference  on  the  dayside,  as  compared  with 
the  nightside  where  photoelectrons  are  absent.  Fig.  4 
shows  the  self-consistent  electric  potential  profiles  for  ^ 
both  polar  wind  solutions.  The  potential  drop  in  the 
dayside  solution  is  about  5  V,  which  is  considerably  j- 

larger  than  that  on  the  nightside  (about  2  V),  and  is 
comparable  to  value  deduced  from  observations  [Win- 
ningham  and  Gurgiolo ,  1982]. 

The  magnitude  of  the  self-consistent  ambipolar  elec¬ 
tric  field  has  a  significant  impact  on  the  ion  outflow 
velocities.  The  0+  and  H+  velocities  in  our  two  po¬ 
lar  wind  solutions  are  shown  in  Fig.  5.  It  is  clear  that 
both  ion  species  have  considerably  smaller  outflow  ve¬ 
locities  in  the  nightside  results,  in  agreement  with  the 
observed  day-night  asymmetry  in  ion  velocities  [Abe 
et  al ,  1993b]. 

The  absence  of  photoelectrons  in  our  nightside  cal- 
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Figure  4.  Self-consistent  ambipolar  electric  potential  pro¬ 
files  in  the  dayside  (solid)  and  nightside  (dashed)  polar  wind 
solutions. 


Figure  5.  Profiles  of  ion  outflow  velocities  in  the  dayside 
(solid)  and  nightside  (dashed)  polar  wind  solutions. 

culations  also  leads  to  the  absence  of  a  noticeable 
electron  anisotropy,  in  contrast  to  the  dayside  situa¬ 
tion.  The  parallel  temperatures  for  the  upward-  and 
downward-moving  electron  populations  in  both  of  our 
polar  outflow  solutions  are  shown  in  Fig.  6.  The  day- 
side  solution  reveals  a  temperature  anisotropy  between 
the  upwardly  and  downwardly  moving  electrons,  *.e. 
Te,up  >  Te}down  (the  ||  subscript  is  omitted  for  simplic¬ 
ity).  This  anisotropy  is  entirely  due  to  the  photoelec¬ 
trons.  In  the  night-time  solution,  the  temperatures  for 
the  two  electron  populations  only  differ  by  about  0.2% 
at  the  altitudes  where  the  satellite  measurements  were 
made  (900  -  1700  km)  [Abe  et  al. ,  1996].  The  day- 
side  and  nightside  electron  temperatures  together  thus 
demonstrate  the  consistency  between  our  photoelectron 
scenario  and  the  day- night  asymmetry  observed  by  the 
Akebono  satellite  [Yau  et  a/.,  1995;  Abe  et  al ,  1996]. 
Because  such  a  temperature  anisotropy  was  observed 
in  the  dayside  polar  wind  [Yau  et  a/.,  1995]  but  seems 
absent  on  the  nightside  [Abe  et  al ,  1996],  the  role  of 


Figure  6.  Parallel  temperatures  for  upwardly  and  down¬ 
wardly  moving  electrons,  the  ||  subscript  has  been  omitted 
for  simplicity.  The  dayside  temperatures  are  represented  by 
the  solid  and  dot-dashed  lines;  the  nightside  temperatures, 
which  are  virtually  equal  for  the  two  electron  populations, 
are  represented  by  the  dashed  line. 

photoelectrons  in  our  model  is  consistent  with  the  ob¬ 
served  polar  wind  scenario. 

4.  Conclusion 

Our  work  on  the  polar  wind  is  motivated  by  the 
increasing  experimental  evidence  that  photoelectrons 
may  affect  the  dynamics  of  the  polar  wind.  In  partic¬ 
ular,  the  day-night  asymmetry  in  the  ion  and  electron 
features  suggests  that  photoelectrons  may  be  the  dom¬ 
inant  effects  in  the  polar  outflow. 

In  order  to  study  the  impact  of  photoelectrons  on  the 
polar  wind,  we  rely  on  a  self-consistent  hybrid  model 
that  provides  a  global  kinetic  collisional  description  for 
the  polar  wind  photoelectrons.  The  model  is  hybrid 
in  that  it  consists  of  a  fluid  part  for  the  thermal  elec¬ 
trons  and  a  kinetic  collisional  part,  and  should  be  dis¬ 
tinguished  from  other  hybrid  schemes  where,  for  exam¬ 
ple,  the  electrons  are  treated  as  a  massless  neutralizing 
fluid.  Specifically,  in  this  model,  photoelectrons  (which 
are  treated  as  test  particles  because  of  their  low  relative 
density),  and  all  the  ion  species  (H+,  0+)  are  described 
by  a  global  kinetic  collisional  approach,  while  thermal 
electron  properties  and  the  ambipolar  electric  field  are 
determined  by  a  fluid  calculation. 

In  order  to  examine  the  anisotropic  kinetic  effects 
of  photoelectrons  as  a  source  of  the  day-night  asym¬ 
metry  in  the  polar  outflow,  we  have  generated  two  po¬ 
lar  wind  solutions,  one  corresponding  to  the  daytime 
outflow,  and  the  other  for  the  nightside.  By  compar¬ 
ing  our  dayside  and  nightside  results,  we  have  demon¬ 
strated  the  consistency  between  our  photoelectron  sce¬ 
nario  and  the  observations.  In  particular,  we  have 
shown  that  the  presence  (absence)  of  photoelectrons  on 
the  dayside  (nightside)  may  explain  the  observed  day- 
night  asymmetry  in  the  polar  outflow.  For  example,  the 
anisotropy  between  the  upward-  and  downward-moving 
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electrons,  which  is  observed  in  the  sunlit  polar  outflow 
but  seems  absent  on  the  nightside,  only  appears  in  our 
day  side  polar  outflow  results.  We  have  also  shown  that 
the  presence  of  photoelectrons  leads  to  a  considerably 
larger  self-consistent  ambipolar  electric  field,  which  cor¬ 
responds  to  a  larger  electric  potential  difference  along 
the  magnetic  field  line.  The  magnitude  of  the  ambipo¬ 
lar  electric  potential  has  a  significant  impact  on  the  ion 
outflow  velocities.  Thus,  both  the  H+  and  0+  veloci¬ 
ties  in  our  results  are  remarkably  higher  on  the  dayside, 
in  agreement  with  satellite  measurements. 
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